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ABSTRACT 

Motivated by the recent investigation of neutrinoless r-lepton decays by the CLEO collab- 
oration, we perform a systematic analysis of such decays in a possible new-physics scenario 
with heavy Dirac/Majorana neutrinos, including heavy-neutrino nondecoupling effects, fi- 
nite quark masses, and quark as well as meson mixings. We find that r-lepton decays into 
an electron or muon and a pseudoscalar or vector meson can have branching ratios close to 
the experimental sensitivity. Numerical estimates show that the predominant decay modes 
of this kind are r~ — > e~0, r~ — > e _ p°, and t~ — ► e~7r°, with branching ratios of order 
1(T 6 . 
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1 Introduction 



Recently, the CLEO collaboration has reported their experimental results on 22 neu- 
trino less decay channels of the r lepton that violate lepton flavour and/or lepton number 
Candidates for lepton-flavour/number- violating events have been found in the decays 
t~ — > e~K*°, t~ — > e~n + K~, t~ — > [i~n~K + , and r~ — > fi + ir~K~. Such decays are 
strictly forbidden in the minimal Standard Model (SM) due to the fact that the light neu- 
trinos, u e , Up, and u T , are exactly massless, so that chirality conservation implies lepton 
number /flavour conservation to all orders of the perturbative expansion. Since there is no 
fundamental theoretical reason for lepton number /flavour conservation in nature, future 
confirmation of the CLEO candidates may point towards physics beyond the SM and, in 
particular, to some modification of the lepton sector. Such possible lepton-number /flavour- 
violating effects can naturally be accounted for in the context of leptoquark models 0, 
left-right-symmetric models 0, imparity- violating supersymmetric scenarios |4[], or theo- 
ries containing heavy Dirac and/or Majorana neutrinos H§ . 

In this paper, we will study the size of new-physics interactions in models with heavy 
Dirac and/or Majorana neutrinos. In such scenarios, decays of the r lepton into three 
charged leptons, such as r — > eee, etc., have been analyzed in Ref. J7J. Here, our main 
interest will be devoted to semileptonic decays of the r lepton. Specifically, we will analyze 
decays of the type t~ — ► e~7r°, t~ — ► e~p°, etc. In a previous work M, three of the 
numerous decay channels of this type were considered in the framework of a theory with 
heavy Dirac neutrinos. We will extend that analysis by including lepton-number-violating 
interactions due to Majorana neutrinos, heavy-neutrino nondecoupling effects |7|i|? finite- 
quark-mass contributions, Cabibbo-Kobayashi-Maskawa (CKM) quark mixings, and meson 
mixings. We will perform a complete analysis, which comprises all the ten decay channels 
of the type t~ — > e~M°, where M° denotes either a pseudoscalar or vector meson. The 
effect of a modified lepton sector on r decays into two mesons, i.e., r" - -> / '.!/, U.; . will 
be estimated in a separate communication |10| . 



In our calculations, we will adopt the conventions and the model described in Ref. || . 
In this minimal model, which is based on the SM gauge group, the neutrino sector is 
extended by the presence of a number (n#) of neutral isosinglets leading to tir heavy 
Majorana neutrinos (Nj), while the quark sector retains the SM structure. If the theory 



contains more than one neutral isosinglet, then the heavy-light neutrino mixing [I I 



3 n R 
1=1 3=1 

scales as {m^ D {rn^) 2 mD)iu where mo is the Dirac mass matrix related to the breaking 
of the SU(2)l gauge symmetry and rriM is a general ur x tir isosinglet mass matrix. 
Light neutrinos (ui) can radiatively acquire masses in compliance with experimental upper 
bounds ||, whereas the lepton-flavour- violating mixings (s^) 2 are dramatically relaxed 
and do not obey the traditional seesaw suppression relation, (s^ ! ) 2 oc rayjm^. Then, the 
(s^) 2 may be viewed as free phenomenological parameters, which may be constrained by 
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a variety of low-energy data fllT| , [T2fl . Throughout this paper, we will consider the following 
conservative upper limits for the flavour- violating mixings [12] : 



(s U l) 2 , (sin 2 < 0.015, (si-) 2 < 0.050, (sl^) 2 (sin 2 < ID" 8 - (1.2) 

Note that these upper bounds are sensitive, to a large extent, to the degree of confidence 
level (CL) considered in the global analyses (i.e., 95% or 99% CL) and further model- 
dependent assumptions 171]. 

In an n^-generation model, the couplings of the charged- and neutral-current inter- 
actions are correspondingly mediated by the mixing matrices 



B ij = *52 v ik u kj, Cij = ^UkiUy, (1.3) 
fe=l fe=l 



where V 1 and U u are the unitary matrices that are needed to diagonalize the charged- 
lepton and neutrino mass matrices, respectively. B and C satisfy a number of identities 
that guarantee the renormalizability of the model [§,0. Such identities are found to be 
very helpful in order to reduce the number of free parameters present in such theories and, 
by the same token, to establish relations between B, C, and the heavy-neutrino masses. 
For definiteness, in our numerical calculations, we will use a model with two right-handed 
neutrinos. In such a scenario, we have |7j 

B lNl = pV4gL ' , B lN2 = ip- 1/A B lNl , (1.4) 



1 + pV2 

where p = m 2 V2 /m 2 Nl , with N\ and N 2 being the heavy Majorana neutrinos. Furthermore, 
the mixings C^jv- are given by 



P n ° v 2 -1/2 

^i-^i = I _|_ pi/2 ' C N2 n 2 = p C NlNl , 

Cn x n 2 — —Cn 2 n 1 = ip CWiATi- (1-5) 

Obviously, our minimal scenario only depends on and m^ 2 — or, equivalently, on 
and p — , and (si) 2 , which are assumed to satisfy the constraints in Eq. ( |1.2| ). 

The outline of this work is as follows. In Section 2, we will calculate analytically the 
branching ratios of the decay processes r~ — > e~M°. Technical details will be relegated to 
the Appendix. Our numerical results will be presented in Section 3. Section 4 contains our 
conclusions. 



2 r" Z'-M° 

Charge conservation forbids the lepton-number- violating decays of a r lepton into a meson 
and an antilepton. For the same reason, the outgoing meson has to be neutral. The recent 
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CLEO experiment |]J observes an event for the decay r _ — > e~K°* within the signal region, 
which is still consistent with the estimated background due to hadron misidentification. The 
same experiment has considerably lowered the upper bounds on the rates of the decays with 
one p° or one K°* in the final state. 

The scattering-matrix element of r~ — > l'~M° receives contributions from 7-exchange 
graphs, Z-boson-exchange graphs, and box graphs, 

S{t~ -> I'-M ) = S^r- -> I'- M°) + 3 z (t- -> l'~M°) + S Box (t- -> l'~M°) . (2.1) 

Feynman diagrams pertinent to these decays are shown in Fig. 1. The 7 and Z-boson ampli- 
tudes factorize into leptonic vertex corrections and hadronic pieces. The loop integrations 
are straightforward. The hadronic parts are local. Exploiting translation invariance, the 
phases that describe the centre-of-mass motion of M° may be isolated and one is left with 
space-time independent hadronic matrix elements. These phases assure four-momentum 
conservation. The 7 and Z-boson amplitudes read 

2 c 2 , „A» 



4M W 

-G:f^^{m T {\ + 75) + m v {\ - 75)) 
' q z 



x(M°|^(0) 7 ^(0) - ^(0) T/1 d(0) - ^(0)7^(0)|0), 
S z ( T -^rM°) = (2n)U^\p-p'-p M )^-Ffu l ,Y(l- l5 )u T 

x({M \u(0)^(l- l5 -^s 2 w )u(0)\0) 

-(M°|d(0) 7M (l- 75-^4^(0)10) 

-<M°|s(0) T/i (l - 75 - ~4>(0)|0>), (2.2) 

where p, p', and j?m are the four-momenta of r, I', and M°, respectively, q = p — p' = pm, 
aw = «em/ sin 2 6w ~ 0.0323 is the weak fine-structure constant, and u(x), d(x), and s(x) 
are quark- field operators acting on the meson states \M°). In Eq. (|2.2| ), , , and G^ 1 
are form factors, which may be found in Ref . [0] . 

The box diagram is more involved, as it contains a bilocal quark operator. Taking 
the difference, X, and the averaged sum of the space-time coordinates of the two hadronic 
vertices as integration variables, using translation invariance, and performing the integra- 
tion over the leptonic space-time coordinates, one arrives at the following expression for 
the box amplitude: 

2 n G +n R 



rv 7r _ r 
S Box (r- ^l'-M°) = (2v)H^\p-p> -p M )^f- £ B ViKi 



(2^ 
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x[ E EK A U/^e- W 

d a ,b=d. ~ ~ ' 1 



x (M°|4(f )7m(1 - 7 5 )4 Ul) P07„(l " 7 5 R( - y) |0) 
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TIG 



Y,v: d y udl J d*x e^ i ~^ x 



x (M°|ti( - |) 7l ,(l - 75 )^)(-X) 7m (1 - 7 >(| ) |0) 
+ (W — > G), (2.3) 

where (W — > G) stands for the terms obtained by replacing one or two bosons with 
unphysical charged Higgs bosons and S F Ut ' dl \x) are the and dj-quark propagators in 
coordinate space. An exact expression for the quark propagator is not known, but using 
the free-quark propagator is known to yield a good approximation for large momentum 
transfers. The corresponding box amplitude which contains only free quarks, receives its 
dominant support from momenta with virtualities comparable to My/ (if the virtual lepton 
is light) or larger. Only a small fraction of the Feynman-parameter space permits poten- 
tially large contributions from low squared-momentum values. Therefore, it is plausible 
to approximate the quark propagator with the free-particle propagator and neglect the 
momenta of the external leptons and the X dependence of the quark wave functions as 
well. Thus, we recover the free-quark expressions for the box functions and evaluate the 
hadronic matrix elements by taking the quark current operators to be local. In this way, 
Eq. (|2.3|) simplifies to 

S Box (r-^l'-M°) = (2n)U^\p-p'-p M )-^u in ,(l- l5 )u T 



Bo r(M°womi -7b)«(o)io> 

-iTl'dadb I 71 rO 



E FroT b <M K(0)7^1-75)4(0)|0>j, (2.4) 

d a ,b=d,s 

where F^ dadh and F^' ^ u may be found in Appendix A. 

To calculate hadronic matrix elements, we invoke the hypothesis of the partial con- 
servation of axial vector currents (PCAC) p^Jl5|Jl6| . 

A p ^x)= i V2f P d,P(x) + ---, (2.5) 

where the dots denote terms not contributing to the meson-vacuum amplitude, fp is the 
decay constant of the pseudoscalar meson P, represented by the field P(x), and A^(x) is 
the axial-vector current having the same quark content as P. The pion decay constant 
is f n = 92 MeV. Furthermore, we exploit the vector-meson dominance (VMD) relation 

2 9 2 

Tfl Til TTlj, 

OH?) = ^p m W + ^ffc^W sin ° v + 2^ 0M(x) cos ° Vl (2,6) 
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and its extension for any vector current 

2 

V?(x) = Ty^W- (2-7) 

In Eq. (|2.6j), jg m (x) is the electromagnetic current, p M (x), and ^(x) are the p-, u;-, 

and 0-meson fields, respectively, j p , and 7^ measure the strengths of their couplings to 
the photon, and 9y is the usual mixing angle of the octet and singlet vector-meson states. 
In Eq. (|2.7|), V^f is the vector field having the same quark content as the vector meson field 



V. Equation Q2.71) is based on the assumption that the dominant contribution to the form 



factors is due to the vector mesons, which works very well for the electromagnetic current 



The calculation of the hadronic matrix elements proceeds as follows. One expresses 
the quark operators, which appear in the hadronic matrix elements, in terms of the axial- 
vector [y4^(x)] and vector currents [V^(a;)] that have the same quark content as the pro- 
duced pseudoscalar (P) and vector mesons (V). Then, one applies Eqs. ( |2.5| ) and ( [2.7|) . 
The relevant matrix elements read 

(Q\A*(x)\M(jp M )) = SMpVtfpppre-**", 

9 

TOP- 

(0\V^(x)\M( PM )) = 5 M y-^-ey^py,\v)e- ip * x , (2.8) 

where e v ^ stands for the polarization vector of the vector boson V, and the Kronecker 
symbols, Smp and S MV , assure that the matrix elements give non-zero contributions only 
if the final-state quantum numbers match those of the vector and axial-vector currents. 
The matrix elements appropriate to mesons in the final state emerge from Eq. ( j2.8|) by 
Hermitean conjugation. 

The decomposition of the vector and axial-vector currents into meson field operators 
depends on the quark content of the meson (for the pseudoscalar mesons, see Table I). 
The quark content of pseudoscalar mesons having zero isospin and zero hypercharge is 
not yet definitely established [|i8| , [T9| . The mixing of SU(3)-octet and SU(3)-singlet meson 
states with zero isospin and zero hypercharge is usually parameterized by some angle, 
9p, which is not precisely known. The corresponding mixing angle for vector mesons is 
called Qy. From the study of (f> decays it is known that 6y is very close to the ideal value 
arctan(l/\/2). Notice that the state \M) and the corresponding field M(x) have opposite 
quantum numbers. This is due to the convention (0\M(x)\M(p)) = e'^ej^ip), i-e., the 
meson field annihilates the corresponding meson state. In the second line of Table I, we 
indicate the relevant creation and annihilation operators that are contained in the meson 
states. Here, b u and d s are the annihilation operators of the quark u and the antiquark s, 
respectively, and b\. and d\ are their creation operators. The quark structure of the vector- 
meson states and fields may be read off from Table I after the replacements K — > K ±0 , 
tt ± -> p ± , 7T° -> p°, K° -> K°\ K° -> K°\ 7/8,1 -> 08,i, V -> 0, v' -> and 6 P -> By. 
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Following the procedure outlined above, we obtain the following expressions for the 
■> l'~M° matrix elements: 



T(t- -> l'-K°) 

t(t~ -> r k°) 

T(r- -> l'-7T°) 
T(r- -> Z'-77) 



T(r" - V-fP) 



T{r- -> r0) 



2Q: 



16M 



16M 



1a 



w 



UlH»(l - 76)«r W + <r + iffl, 



2Q: 



IV 



16M 



IV 



,2c, 



+ 



/ C P y/2s P 



KBox + {—7= 



V2s P 



Box 



>/3 + >/3 



2Q: 



16M W 



2 
11 

%/2c P 



. V2c P 2s 



V3 V3 



)*z 



V2cp _ 2sp Vss 
1 V3 V3> Box 



S P s W rl'uu , / V^CP Sp VM 

+ - =: '7,w + I ^ + ^ ]/ 7;.u- 



2tt 



IV 



uv^^l - 75) 



1/,, 



M rprl'sd 



TiTl'ds 



■ 2 2 
l0i W - /1 \ m K°* u r i 



4 m PQ ^ 



16M^ 7po ^ 



F; V (7, -%/)(! -75) 



_ G ?'^(m T (l + 75 )+m,(l- 75 
+%7m(1 - 7sK [awF* + i<'r + jiffl}, 



16M 2 



iv 70 



<(7m-^)(1"75) 



_ G ^'!^?l( mT(1 + 75) + ro| , (1 _ 75 

+ ^(i- 75 K[(^ + ^)^ + ( ^-^)<r 



4 m 5 cM 



{ 



{ 2Vs VE } Box { Vs VQ 



1s 2 w s v _ 



f; 1 (7m 



% 4 



)(l-75) 



7 



-G 



m T (l + 75) + m v (l -75) 



+iii'7 M (l - 7 5 )m t [(-^=c 2W ^JK 



/ S y 1 CV ^ Tprl'uu 



2V3 Vq' Box { Vs VE } Box ^'' 



[2.9) 



where we have introduced the short-hand notations sw = sin^iy, cw = cos^jy, c 2 w — 
cos2#n/, and similarly for 9p and Qy. 

The branching ratios for pseudoscalar mesons can be compactly written in the form 



B(t~ -> l'~M°) 



1 m T \2(ml, mf,,m 2 M0 ) 2 (m 2 - m 2 ,) 2 - m 2 m {m 2 T + m 2 ) 



|a M o 



(2.10) 



87r T T m 2 

where the form factors a M o are listed in Appendix A, Y T = 2.16 10~ 12 GeV is the total 
width of the r lepton measured experimentally, and X(x, y, z) = (x—y — z) 2 —Ayz. Similarly, 
in the case of the vector mesons, one finds 



1 m T A2 (m 2 , mf,, m 2 M o) 



8tt r T 



mi 



\CM°\ 



>| 12m 2 , 



m 



(ml + mf,) (2(m 2 T - mf,) 2 - m 2 MQ (m 2 T + mf,) - m 4 M0 ) 



m 2 m^ 

,2\2 



} 



(m 2 - m 2 ) 2 + m 2 M0 (m 2 T + mf,) - 2m\ l0 



2 > >o 2 



mtm 



M" 



-63te{(a M ° + b M o)c* M a} 



(m 2 - mf,) 2 - m 2 M0 (m 2 + mf) 



m 2 T m 2 M0 



(2.11) 



where a M o, b M o, and c M o may also be found in Appendix A. 



3 Numerical results 

In our numerical analysis, we will assume that the SM is extended by two right-handed 
neutrinos, as described in the Introduction. The additional parameters in this scenario 
are the two heavy-neutrino masses, m^ 1 and m^ 21 and the three mixing angles, s v £, s v £ , 
and s v £ . These are free parameters of the model, which may be limited by experiment. 
The upper bounds on s v £, s v £ , and s v £ are given in Eq. ( |1.2|) . On the other hand, the 
perturbative unitarity relations, 

' ' < ~ (3.1) 



m Ni 2 

lead to a global upper bound on m^, 
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where p is defined after Eq. (|1.4|) and it is understood that p > 1. In this context, we 



should mention that, adapting the results of Ref. |20| based on a renormalization-group 
analysis in a four-generation Majoron model, one may find a bound which is slightly more 
restrictive than Eq. ( |3.2|) but still lies in the same ball park. 

Furthermore, our results depend on hadronic observables and quark-level parameters 
such as the CKM-matrix elements, the quark and meson masses, the mixing angles of the 
meson singlet and octet states, the pseudoscalar-meson decay constants, and the coupling 
strengths of the vector mesons to the gauge bosons. In our calculations, we use the max- 



imum experimental values for the CKM-matrix elements |L8| and the quark-mass values 

00 



m u = 0.005 GeV, m d = 0.010 GeV, m s = 0.199 GeV, 

m c = 1.35 GeV, m b = 4.3 GeV, m t = 176 GeV. (3.3) 

We keep all quark masses finite, since, e.g., the c-quark and t-quark contributions to the 
box amplitudes turn out to be comparable. The mixing angle for vector-meson nonet states 
may be determined from the quadratic Gell-Mann-Okubo mass formula to be 0y = 41.3°. 
We treat Op as a free parameter because its value is not yet well established [ 19|,p4| . For 
the most part, we use Op = —23°, the value extracted from e + e~ — > e + e~77* — > e + e~(P — > 
77) experiments f22|. This value is consistent with a previous analysis [19|]. For the 



pseudoscalar-meson decay constants, we use the experimental values fl8|,|2~ 



fa± = 92.4 MeV, f K ± = 113 MeV, 

fa = 84.1 MeV, /„ = 94 MeV, fa = 89.1 MeV, (3.4) 
and exploit SU(3) flavour symmetry, 

fa° = fn° « /k±- (3.5) 

The constants jy are partly extracted from the V — > e + e~ decay rates, with the result that 

7p o = 2.519, j u = 2.841, l4> = 3.037, (3.6) 

and partly estimated assuming SU(3) symmetry: we put 7^.0 = 7 p o because K*° and p° 
are members of the same SU(3) octet, while and u are mixtures of octet and singlet 
states. Notice that all 7^ values in Eq. ( |3.6|) are very similar in size. 

Having specified our input parameters, we will now discuss our numerical results. 
The widths for the decays with K°, K°, K *, or K°* in the final state only receive contri- 
butions from box diagrams. The branching ratios for these decays are found to be always 
smaller than 10~ 14 , that is, much smaller than present experimental sensitivities (~ 10 -6 ), 
rendering these decay modes uninteresting from the experimental point of view. Therefore, 
we will not pursue their study any further. 

For definiteness, we will consider decays of the form r~ — > e~M° only — we set (s^ 1 ) 2 ~ 
to satisfy the third inequality in Eq. ( |1.2|) . Of course, our estimates are also valid for 
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the t~ — > pT M° decays with {s u £) 2 = 0. Our results for the branching ratios B(t~ — > 
e^7i° /rj/r]' / p°/<f)/uj) are illustrated in Figs. 2-6. Each figure describes the dependence of the 
branching ratios on two of the free parameters, one varied continuously and the other one 
in a discrete manner. All other parameters are kept fixed. Figure 2 shows the dependence 
of the branching ratios on itln = itln 1 = rriN 2 and (s u £) 2 . The most promising modes are 
t~ — > e~<f), r~ — > e~p°, and t~ — > e~7r°, which, for maximum values of mjy and (s v £) 2 , 
reach branching fractions 



B(t~ 


— > e 0) ~ 


1.6 


10 


B(t~ - 




0.9 


10 


B(t~ - 




1.0 


10 



(3.7) 

This has to be compared with the present experimental bounds [ftl,^ 

B(r~ -> e~p°) < 4.2 • 10~ 6 , 
^(r-^e-Tr ) < 1.4 • 10" 4 , 

B(t~ -> ^7T°) < 4.4 -10" 5 , (3.8) 

at the 90 % CL. Unfortunately, B{t~ —>■ e~<p) has not been measured yet. We conclude that 
an experimental investigation of t~ — ► e _ and a more precise determination of .B(t _ — » 
e~7r°) and B(t~ — ► /U~7r°) would be highly desirable. In the high-m^v! limit, Fig. 2 shows 
the quadratic m^v dependence for all branching ratios, except for B(t~ — > e^cj). In the 
't Hooft-Feynman gauge, this behaviour originates mainly from the Z-boson amplitudes, 

?z ■ 

At this stage, some important comments are in order. For fixed {s u £) 2 values, N% 
and N2 do not decouple from our theory as their masses become large as compared to 
Mw [[OH- As has been mentioned in the Introduction, s v £ oc mfl/m« oc mn/m^^ and this 
nondecoupling feature can be traced to the large SU(2)l Dirac components, mo, present 
in our model ||. Obviously, if we fix mp and take the limit — > 00, the heavy neutrinos 
will decouple from our low-energy processes, leading to vanishing effects This will be 



illustrated in greater detail in Figs. 4 and 5. However, for heavy neutrinos, with masses 
in the 1-10 TeV range, there will be an interesting nondecoupling "window" arising from 
potentially large Dirac mass terms mjj. It is precisely this nondecoupling "window" which 
we are exploiting here to make our effects sizeable. 

In the case of B(t~~ — > e^ir jr\j K°), we recover the expressions of Ref. || for the 
Dirac-neutrino scenario if we omit the nondecoupling terms proportional to m 2 N in Eq. ( |A.2| ) 
and in the Z-boson-mediated amplitudes. The results of Ref. || are comparable to ours for 
m N of order M w , but they fall short of our results by up to factor of 50 for m N in the TeV 
region. In the case of uo production, there is a destructive effect between logarithmic and 
quadratic m^-dependent nondecoupling terms coming from photon and Z-boson-mediated 
amplitudes, respectively, and meson-mixing effects, which show up as a minimum of the 
branching ratio for ~ 1.6 TeV. 

We now turn to genuine Majorana-neutrino quantum effects. Figure 3 displays the de- 
pendence of the branching fractions on the ratio m^/m^ for the fixed values = 1 TeV 
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and 0.5 TeV. We emphasize that, just like in the lepton case |7|], B{t~ — > er n° / r\ / r( / p° / '(f>) 
assume their maximum values for m^/m^ = 2-4 rather than in the Dirac scenario, 
m = mN 2 - The only exception is the decay r~ — > e~uj, where the maximum value 
is shifted to larger values of m^/^Vi, of order 20, due to the accidental cancellations 
mentioned above. 

Figures 4 and 5 illustrate the dependence of the branching ratios on {s v £) 2 and (s£ e ) 2 
in the heavy-Dirac- neutrino scenario with = mjv 2 =4 TeV. As may be seen in Fig. 4, 
the {s v £) 2 dependence of B{r~ — > e~ n° / r] / 'rf / p° / <p / 'uj) is quadratic over the most part of the 
{s v £) 2 range and for any of the (s^ ) 2 values considered. From Fig. 5 we see that the (s^ ) 2 
dependence of B{t~ — > e~ n° / r] / 'rf / p° / '</> / 'uj) is approximately linear for (s u £) 2 < (s^) 2 , 
while it becomes quadratic for (s^ ) 2 > (s^j) 2 . The (s^ T ) 2 and (s£ e ) 2 dependences studied 
above are closely related to the decoupling behaviour of the isosinglet scale uim- As we have 
emphasized above, in the limit —>■ oo for constant — or, equivalently, for constant 
mjy and vanishing tud, i-e., for (s^) 2 — > — we should recover the decoupling limit, where 
the branching ratios vanish as the isosinglet mass terms mu are sent to infinity |24j]. It is 



then evident that the aforementioned (non) decoupling "window" is directly related to the 
SU(2)l Dirac terms mo and is reflected in the actual {s v £) 2 and {s 1 ^) 2 dependences seen 
in Figs. 4 and 5. 

In Fig. 6, we plot B(t~ — > e~r]/r]') versus = TOjvi = ^n 2 i assuming in turn the 
unmixed case {6p = 0) and Op = —10°. We see that, for Op decreasing, B{t~ — ► eTr() 
increases considerably, while B{r~ —>■ e~i]) grows just slightly. This illustrates that it is 
important to allow for nonvanishing Op in realistic calculations. It is also interesting to 
observe that, if tanOy = 1/(v / 2c2vk), the dominant nondecoupling terms proportional to 
m 2 N are quenched in B{t~ — > 1~uj). 



4 Conclusions 

Motivated by the recent experimental search for lepton-number /flavour- violating semilep- 
tonic r-lepton decays which are strictly prohibited in the SM, we have explored the 
potential of extensions of the SM by heavy Dirac and/or Majorana neutrinos to account for 
t~ — > l'~M° decays, where I' = e, p and M° is a neutral pseudoscalar or vector meson, with 
branching ratios which are in line with the experimental results. Since such models predict 
appreciable branching fractions for lepton-flavour /number- violating leptonic decays of the 
t lepton 0, they are also expected to be promising candidates for explaining the problem 
at hand. In fact, we have found branching fractions in excess of 10~ 6 for the channels 
t~ — > e~(p, t~ — ► e~p°, and r~ — > e~7T°. Our value for B{r~ — ► e~p) is comparable to 
the present experimental sensitivity Unfortunately, the experimental upper limit on 
B{t~ — > e _ 7r°) still exceeds our result by two orders of magnitude ||23|| . For some reason, 
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the decays r _ — > or r _ — > which prevail in our numerical estimates, have not yet 
been studied experimentally. At this point, we should like to encourage our experimental 
colleagues to undertake a search for this decay channel. 

An important feature of our model is that the t~ — > l'~M° decay amplitudes exhibit 
a quadratic dependence on the heavy-neutrino masses, and m^. This nondecoupling 
dependence is closely related to the large SU(2)i-breaking Dirac terms that are allowed 
to be present in our minimal t/iree-generation seesaw-type scenario ]7||J. These m? N . terms 
are negligible for neutrino masses below 200 GeV, but they are dominant in the TeV region, 
where they may lead to an enhancement by a factor of 50 of the respective analysis with 
these terms omitted [|J. The same nondecoupling terms give rise to a m N2 /m Nl dependence 
of the r~ — > l'~M° decay amplitudes which is similar to the one encountered for the decays 
t~ — > e + e~e~ , etc. [[7]. In particular, semileptonic branching ratios take their maximum 
values for m^/m^ ~ 2-4. The r _ — > 1'~oj decay rate is unobservably small due to 
a destructive meson-mixing effect, which considerably screens the dominant Z-exchange 
interaction. 

The extension of the vector-meson dominance hypothesis to general vector currents 
has enabled us to calculate the decays with vector mesons in the final state. The quark 
content of meson wave functions, which, for instance, is reflected in the mixing angles, 
6p and 9y, is also important. We have illustrated this for the production of r], rf, and (ft 
mesons. 
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A Appendix 



The amplitudes a^o, b^o, and ca/o appearing in Eqs. (|2.10|) and ( |2.11|) may be decomposed 
into the form factors F%' * ad \ F^ u , Ff' , F^ 1 ' , and G*' in the following way: 
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(A.1) 



The form factors F Box adh , F B ^ X U , Fp' , F^', and G^', which also appear explicitly in Sec- 
tion 2, may in turn be decomposed into elementary vertex and box functions, F~, G 7 , F z , 



Fboxi an d R Box- The form factors Fp , F^ , and G^ together with the elementary loop 
functions F 7 , G 7 , F Z} and F Box may be found in Ref. 0. Here, we list F Box adb and F B ' , 
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To our knowledge, the box function i?B oa; may not be found elsewhere in the literature. 
After a straightforward calculation, we obtain 
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For the reader's convenience, we evaluate H Box for some special arguments 
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Figure Captions 



Fig. 1: Feynman graphs pertinent to the semileptonic lepton-flavour- violating decays 
T~ -> l'-M°. 

Fig. 2: Branching ratios versus heavy-neutrino mass = = m^ 2 for the 

decays r~ — > e~7r° (solid line), t~ — > e^r] (dashed line), r~ — > e - ?/ (dotted 
line), t~ — > e~p° (dot-dashed line), r~ — > e~0 (dot-long-dashed line), and 
T~ -> e~cj (thick-dotted line), assuming (s£) 2 = 0.01, 0.02 < (s£ T ) 2 < 0.05, 
9 P = -23°, and y = 41.3. 

Fig. 3: Branching ratios versus ratio m N , 2 /m Nl for the decays of Fig. 2, assuming 
m Nl = 1 TeV (0.5 TeV), (s£) 2 = 0.01, (s£) 2 = 0.05, # P = -23°, and 
9 V = 41.3. 

Fig. 4: Branching ratios versus {s v £) 2 for the decays of Fig. 2, assuming rajvi = 
m N2 = 4 TeV, 0.005 < (s£) 2 < 0.014, # P = -23°, and 6 V = 41.3. 

Fig. 5: Branching ratios versus (s u £) 2 for the decays of Fig. 2, assuming = 
m N2 = 4 TeV, 0.01 < (s^f < 0.04, 6 P = -23°, and 6 V = 41.3. 

Fig. 6: B(t~ — > e — 77/77') versus mjy = mjvi = ^a^; assuming (s£ e ) 2 = 0.01 and 
(s^) 2 = 0.04: B(t~ -> e-77) (solid line) and £(7-- -> e"?/) (dashed line) 
for the unmixed case 9 P = 0; £?(t~ — > e~?7) (dotted line) and £?(t~ — > e - ^') 
(dot-dashed line) for the mixed case with 9 P = —10°. 
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Table I: Quark content of the pseudoscalar meson states \M) and field operators M(x). 
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